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Abstract 

Covariant density functional theory, which has so far been applied only within the framework 
of static and time dependent mean field theory is extended to include Particle- Vibration Coupling 
(PVC) in a consistent way. Starting from a conventional energy functional we calculate the low- 
lying collective vibrations in Relativistic Random Phase Approximation (RRPA) and construct an 
energy dependent self-energy for the Dyson equation. The resulting Bethe-Salpeter equation in 
the particle-hole (ph) channel is solved in the Time Blocking Approximation (TBA). No additional 
parameters are used and double counting is avoided by a proper subtraction method. The same 
energy functional, i.e. the same set of coupling constants, generates the Dirac-Hartree single- 
particle spectrum, the static part of the residual p/i-interaction and the particle-phonon coupling 
vertices. Therefore a fully consistent description of nuclear excited states is developed. This method 
is applied for an investigation of damping phenomena in the spherical nuclei with closed shells 
208 Pb and 132 Sn. Since the phonon coupling terms enrich the RRPA spectrum with a multitude of 
p/i£3phonon components a noticeable fragmentation of the giant resonances is found, which is in full 
agreement with experimental data and with results of the semi-phenomenological non-relativistic 
approach. 

PACS numbers: 21.10.-k, 21.60.-n, 24.10.Cn, 21.30.Fe, 21.60.Jz, 24.30.Gz 
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I. INTRODUCTION 



The theoretical analysis of nuclear dynamics by microscopic methods remains undoubtedly 
a crucial issue of present-day nuclear physics. Besides the merely theoretical interest it 
is stimulated by new experimental opportunities for the study of the nuclear many-body 
system. Along with new possibilities to investigate more precisely nuclear systems being 
considered well known, the recent generation of radioactive-beam facilities enables us to 
examine exotic systems with extreme spin and isospin values. For a successful description 
of nuclei close to the limits of stability a powerful and reliable theory should be provided. It 
should be based on a consistent treatment of both ground and excited states and it should 
allow for predictions of nuclear properties in areas, which are hard or impossible to access 
by future experiments. 

Ab initio calculations and multi-configuration mixing within the shell-model can, so far, only 
be applied in light nuclei. At present, for a universal description of nuclei all over the periodic 
table, Density Functional Theory (DFT) based on the mean-field concept provides a very 
reasonable concept. DFT has been introduced in the sixties in atomic and molecular physics 
l|, |2| and shortly after that in nuclear physics, where it has been called density dependent 
Hartree-Fock theory 3, J]. Today it is widely used for all kinds of quantum mechanical many- 
body systems. Density functional theory can, in principle, provide an exact description of 
the many-body dynamics, if the exact density functional is known, but for such systems as 
nuclei one is far from a microscopic derivation and one is forced to determine the functional 
in a phenomenological way. Starting from basic symmetries the parameters are adjusted to 
characteristic experimental data in finite nuclei and nuclear matter. 

One of the most successful schemes of this type is Covariant Density Functional Theory 
(CDFT). It is based on Lorentz invariance connecting in a consistent way the spin and 
spatial degrees of freedom in the nucleus. Therefore, it needs only a relatively small number 
of parameters which are adjusted to reproduce a set of bulk properties of spherical closed- 
shell nuclei (5, 6]. A large variety of nuclear phenomena have been described over the years 
within this kind of model: ground state properties of finite spherical and deformed nuclei 
all over the periodic table J7] from light nuclei Q| to super-heavy elements jol, Q, from 



the neutron drip line, where halo phenomena are observec 



111 ], to the proton drip line [1 



with nuclei unstable against the emission of protons [131]. Rotational bands are treated 
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within the relativistic cranking approximation and the Relativistic Random Phase 

Approximation (RRPA) [16| and quasi-particle RRPA 17] have been formulated as the small 
amplitude limit of time-dependent Relativistic Mean-Field (RMF) models for a description 
of collective and non-collective excitations. This method is successful, in particular, for the 
description of the positions of giant resonances and spin/isospin-excitations as the Isobaric 
Analog Resonance (IAR) or the Gamow- Teller Resonance (GTRj_18|. Recently it has 
been also used for a theoretical interpretation of low-lying dipole 17] and quadrupole 19| 
excitations. 

With a few exceptions, where the quadrupole motion has been studied within the relativistic 
Generator Coordinate Method (GCM) 201 ]. applications of covariant density functional the- 
ory to the description of excited states are limited to relativistic RPA, i.e. to configurations 
of lpl/i-nature. None of these methods, however, can be applied to the investigation of 
the coupling to more complicated configurations, as it occurs for instance in the damping 
phenomena causing the width of giant resonances. It is also known, that in such theories, be- 
cause of the low values of the effective mass, the level density of the single-particle spectrum 
in the vicinity of the Fermi surface is considerably reduced as compared to the experimental 
data. 



Based on the Fermi Liquid Theory (FLT) of Landau [2l| the Theory of Finite Fermi Systems 
(TFFS) of Migdal 22] is another very successful method for the description of low-lying 
nuclear excitations 23], which is known since the fifties. It has several general properties in 
common with density functional theory. First, both theories are know to be exact, at least 
in principle, but in practice the parameters entering these theories have to be determined 
in a phenomenological way by adjustment to experimental data. Second, both theories are 
based in some sense on a single-particle concept. Density functional theory uses the mean 
field concept with Slater determinants in an effective single-particle potential as a vehicle to 
introduce shell-effects in an exact theory and Fermi liquid theory is based on the concept of 
quasi-particles obeying a Dyson equation, which are defined as the basic excitations of the 
neighboring system with odd particle number. Third, in practical applications both theories 
describe in the simplest versions the nuclear excitations in the RPA approximation, i.e. by 
a linear combination of p/i-configurations in an average nuclear potential. 
However, there are also essential differences between these two concepts. First, in contrast 
to DFT, TFFS does not attempt to calculate the ground state properties of the many-body 
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system, but it describes the nuclear excitations in terms of Landau quasi-particles and their 
interaction. Therefore the experimental data used to fix the phenomenological parameters 
of the theories are bulk properties of the ground state in the case of density functional 
theory, and properties of single-particle excitations and of the collective excitations in the 
case of finite Fermi systems theory. Second, in DFT the mean field is determined in a self- 
consistent way and therefore the RPA spectrum contains Goldstone modes at zero energy. 
This is usually not the case in TFFS calculations, which are based on a non-relativistic 
shell-model potential, whose parameters are adjusted to the experimental single-particle 
spectra. Therefore there is no self-consistency in the RPA calculations of TFFS and the 
Goldstone modes do not separate from the other modes. They are distributed among the 
low-lying excitations. Third, modern versions of TFFS go much beyond the mean field 
approximation. Using Green's function techniques the coupling between quasi-particles and 
phonons has been investigated. Based on the phonons calculated in the framework of the 
RPA one has included particle-phonon coupling vertices and an energy dependence of the 
self-energy in the Dyson equation 



24 



251 ]. This leads also to an induced interaction in 



the Bethe-Salpeter equation caused by the exchange of phonons which also depends on the 
energy. The coupling of particles and phonons has also been derived from Nuclear Field 



Theory (NFT) and its extensions 



26 



27, 



. Many aspects of the coupling between the 



quasi-particles and the collective vibrations have been investigated with these techniques 



29 
39 



30 



31 



32 



33 



34, 



35 



36 



37 



38) as well as with other kinds of approaches beyond RPA 



401 ] over the years. 

The present work makes the attempt to find a combination of the basic ideas of density 
unctional theory and Landau-Migdal theory. The concept is similar to earlier work in Refs. 



41 



42 



431 ] . where specific non-relativistic energy functionals have been used to construct 



Self- Consistent Theory of Finite Fermi Systems. We start here from a covariant density 
functional E[p] widely used in the literature. It is adjusted to ground state properties of 
characteristic nuclei and, without any additional parameter, it provides the necessary input 
of finite Fermi systems theory, such as the mean field and the single-particle spectrum as well 
as an effective interaction between the p/i-configurations in terms of the second derivative 
of the same energy E[p] with respect to the density. Thus the phenomenological input 
of Landau-Migdal theory is replaced by the results of density functional theory. The same 
interaction is used to calculate the vertices for particle- vibration coupling 44] . We then apply 
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the techniques of Landau-Migdal theory and its modern extensions to describe vibrational 
coupling and complex configurations. In this way we obtain a fully consistent description of 
the many-body dynamics. 

Two essential approximations are used in this context: First, the Time- Blocking Approxi- 
mation (TBA) blocks in a special time-projection technique the lpl/i-propagation through 
states which have a more complex structure than lpl/i®phonon. The nuclear response can 
then be explicitly calculated on the lplh + lpl/i®phonon level by summation of an infinite 
series of Feynman's diagrams. This approximation has been introduced in Ref. 34j. Ini- 
tially it was called Method of Chronolo gica l Decoupling of Diagrams (MCDD) and it is also 



discussed in recent review articles 



35 



361 ] . Second, a special subtraction technique guar- 



antees, that there is no double counting between the additional correlations introduced by 
particle-vibration coupling and the ground state correlations already taken into account in 
the phenomenological density functional. These two tools are essential for the success of the 
present method. TBA introduces a consistent truncation scheme into the Bethe-Salpeter 
equation and without it it would be hard to solve the equations explicitly. The subtraction 
method is the essential tool to connect density functional theory so far used only on the level 
of mean field theory, i.e. on the RPA level, with the extended Landau-Migdal theory, where 
complex configurations are included through particle-vibration coupling. 
As discussed above there are several versions of density functional theory in nuclear system, 
n the present work we use covariant density functional theory with the parameter set NL3 



451 ] . A first attempt in this direction has been carried out in Ref. [44| , where we treated the 



level density problem in the framework of covariant particle-vibration coupling. 
In Section II we discuss shortly the general formalism of covariant density functional theory, 
we introduce the concept of the energy dependent self-energy S(e) and the vertices of particle- 
vibration coupling in the relativistic framework, and we discuss the response formalism and 
the time blocking approximation for the response function. In Section III we show, how to 
solve the resulting integral equations in detail and present recent numerical applications for 
the spreading width of the doubly magic nuclei 208 Pb and 132 Sn and in Section IV we include 
a brief summary and an outlook for future applications. Appendices A-D provide some more 
details for the formalism used in the calculations. 
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II. FORMALISM 



A. The relativistic mean-field approach 



- d v pn , 

- d u A». (2) 



In relativistic mean-field theory the nucleus is considered as a system of nucleons interacting 
by the exchange of effective mesons, which characterize the quantum numbers of the various 
fields. Conventional RMF theory uses er, uj and p- mesons and the photon as a minimal set 
of bosons which is necessary to generate realistic nuclear fields [5|, |46|, |47j. This concept is 
expressed by the following Lagrangian density depending on the nucleon spinor ip as well as 
on the mesons a, u^, p^ and the electromagnetic fields A 11 : 

£ = V^ 7 ^ t - m)i) + l -d»ad p a - ^m^ 2 



- -ara^ + -m 2 ^ - -htr^ + -m 2 p p»p p - -f^f, u 

- i> (r CT a + r> M + f/% + T^A^j i>. (1) 
Here m is the bare nucleon mass, m a , and m p are the corresponding meson masses, and 



^nv, R^v, F pu are field tensors: 



Fnv duAy 



As usual arrows denote isovectors and boldface symbols will be used for vectors in ordinary 
space. The vertices V entering the interaction term of Eq. ([1]) read 

T a = g a , r£ = g u >f, f/ = g p f^, T» = e^^f, (3) 

where g a , g^, g p and e are the corresponding coupling constants. For the sake of simplicity, 
we use in the following discussions the linear version ([1]) of the Lagrangian. However, it is 
well known that one needs an additional density dependence in order to obtain a reliable 
description of nuclear surface properties. It is either expressed by a non-linear self-interaction 
between the mesons or by an explicit density dependence of the coupling constants g a , g U) 
and g p . In all the numerical applications we use such modifications. They are discussed in 
detail in Appendix [A] 
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The classical variation principle being applied to the Lagrangian density ([I]) leads to the 
following system of coupled equations of motion for the Fermions 

(17^ ~m-J2 T m m (r, t))Vi(r, t) = (4) 

m 

and the mesons: 

(□ + m^)0 m (r, t)) = T J2 ^( r > ^) r m^(r, *). (5) 

i 

Here the index i numerates the nucleons and the index m = {cr, cu, p, e} runs over the set 
of meson fields <p m = {a, u 11 , p^, A 11 }, the corresponding masses and T m are the 

vertices ([3]). The —sign in Eq. ([5]) holds for scalar fields and the +sign for vector fields. 
The set of Eqs. (|4]) and ([5]) defines the standard RMF model. It implies the following four 
approximations: 

i) only the motion of the nucleons is quantized, the meson and electromagnetic degrees 
of freedom are described by classical fields. The nucleons move independently in these 
classical fields, i.e. residual two-body correlations are disregarded and the many-nucleon 
wave function is a Slater determinant at all times. 

ii) in time-dependent applications we use the instantaneous approximation, i.e. the time 
derivatives in the Klein-Gordon equations and consequently the retardation effects for 
the meson fields are neglected although the Dirac spinors ipi as well as the fields m are 
functions of coordinates and time. This is justified by the relatively large meson masses and 
the corresponding short range of the meson exchange interaction. 

iii) Fock terms are neglected. Because of the short range character of the meson exchange 
forces they behave similar to zero range forces, where the Fock terms have the same form as 
the direct terms and contribute only in a change of the relevant strength parameters. Since 
our strength parameters are adjusted to experiment, Fock terms for the heavy mesons are 
implicitly taken into account to a large extend. 

iv) the no-sea approximation means, that nucleon states in the Dirac sea with negative 
energies do not contribute to the densities and currents: the sum over i in the Eq. ([5]) 
includes only occupied levels with positive energy in the Fermi sea. This means we do not 
consider vacuum polarization explicitly. Since our parameters are adjusted to experimental 
data such effects are taken into account in an implicit way. 
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B. The energy-dependent nucleon self-energy 

Let us now consider in detail the Dyson equation, which describes the motion of a single 
nucleon in the nuclear interior: 

(17"^ - m - E)^ = 0. (6) 

Here the total self-energy E contains all the forces acting on the nucleon. As long as one stays 
within RMF approach the nucleon self-energy contains only static and local contributions 
from the mesons and from the electromagnetic fields. 

When we go beyond the mean-field approximation, we have to take into account that in 
the general case the full self-energy E is non-local in the space coordinates as well as in 
time. This non-locality means that its Fourier transform has both momentum and energy 
dependence. Let us now decompose the total self-energy matrix into two components, a 
static local and an energy dependent non-local term: 

E(r, r'; e) = E(r)<5(r - r') + E e (r, r'; e), (7) 

where 

E(r) = ^r m m (r) (8) 

m 

is the static part known from RMF theory in Eq. (j4]). The upper index "e" in E e indicates 
the energy dependence. The Dyson equation for single-particle Green's function reads: 

(e-h v - (3E e {e))G{e) = (3, (9) 

where h v denotes the Dirac hamiltonian in Eq. (j3J) with the energy-independent mean field: 

hP = ap + [3(m + t s ) + E , a = £7, (3 = 7 . (10) 

Most of the static calculations are done for cases, where time-reversal symmetry is valid. In 
all these cases there are no currents in the nucleus and, thus, only scalar E s and time-like 
So components of E survive in the static Dirac hamiltonian of Eq. ( TlOl . In the shell-model 
Dirac basis {1^)} which diagonalizes the energy-independent part of the Dirac equation 

= e k \^ k ) (11) 
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one can rewrite Eq. as follows: 

J2i( £ - £ ^ - Z e kl (e)}G lk ,(e) = 6 kk ,, (12) 
i 

where 

= J d 3 rd 3 r' ^ fc (r)S e (r, r'; e)^(r') , (13) 
G fe ,00 = / dVdV ^(r) /3G(r, r'; £ )/?^(0 • (14) 



The letter indices k, k', I denote full sets of the single-particle quantum numbers. In the 
present work we consider spherical symmetry, so that the spinor \ip k ) is characterized by the 
set k = {(k),mk} and (k) = {n k , jk,^k,^k} with the radial quantum number n k , angular 
momentum quantum numbers jki^k, parity ii k and isospin r k : 

If A ( f(k)(r)^i kjk mM^ s ) | f A nK , 
M r ,s,t) = XT k (t), (15) 

where s and t are the coordinates for spin and isospin. The orbital angular momenta l k and 
l k of the large and small components are determined by the parity of the state k: 



h = jk + |, h= jk - \ for vr fc = (-l) Jfc+ 2 

i (1°) 
h = jk - |, k= jk + \ for n k = (-l) Jfc -2, 

f(k)( r ) an d 9(k)( r ) are radial wave functions and $y TO (#, <p, s) is a two-dimensional spinor: 
$z jm ($, tp, s) = ^(~m s lm l \jm)Y lmi {'&,v)Xrn s {s)- (17) 
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As in Ref. 



441 ] we use the particle-phonon coupling model for the energy dependent part of 
the self-energy S e (e). In the Dirac basis of the spinors 1^) defined in Eq. (fTTj) it has the 
form 

fi,n ' v 

where r\ — ■* +0. In Fig. [T]we show a graphical representation for £^(e). The index n runs 
over all single-particle states in this basis. We distinguish occupied states in the Fermi sea 
\h) (hole states) with cr^ = —1, unoccupied states above the Fermi level \p) (particle states) 
with Op = +1, and states in the Dirac sea \a) (anti-particle states) with negative energies. 
Due to the no-sea approximation the orbits with negative energies \a) are left unoccupied 
and therefore we have a a = +1. However, as it has been shown in Ref. the numerical 

10 




FIG. 1: The particle £y p/ , and the hole ^%i h n components of the relativistic mass operator in 
the graphical representation. Solid and dashed lines with arrows denote one-body propagators for 
particle (p), hole (h), and antiparticle (a) states. Wavy lines denote phonon (//) propagators, empty 
circles are the particle-phonon coupling amplitudes 7^. See text for the detailed explanation. 

contribution of the diagrams with intermediate states a with negative energy is very small 
due to the large energy denominators in the corresponding terms of the self-energy (lisp . The 
index [i in Eq. ( 1181) labels the various phonons taken into account. f2 M is their frequency, 
p M is their transition density and the phonon vertices 7^ determine their coupling to the 
nucleons: 

k'l' 

Vki',ik' denotes the relativistic matrix element of the residual interaction. It is a obtained as 
functional derivative of the relativistic mean field with respect to relativistic density matrix 
P 

V kV ,w = -y"^. (20) 

Opik 

We use the linearized version of the model which assumes that p M is not influenced by the 
particle-phonon coupling and can be computed within the relativistic RPA approximation. 
As long as we neglect the non-linear self-couplings of the mesons or the density dependence 
of the coupling constants g m in Eq. ([3]) the matrix elements V kV lk , in Eq. f fT9|) have simple 
Yukawa form. However, such effects have to be taken into account to guarantee full self- 
consistency of the RRPA equations and we find more complicated expressions, which are 
discussed in detail in Appendix \K\ In the applications investigated in section III we use the 
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non-linear parameter set NL3 45j for the solution of the static Dirac equation ffTTl) as well 
as for the calculation of the matrix elements (120p . 



C. The response function 

Nuclear dynamics of an even-even nucleus in a weak external field is described by the linear 
response function i?(14, 23) which is the solution of the Bethe-Salpeter equation (BSE) in 
the ph channel. In the beginning it is convenient to consider this equation in the time 
representation. Let us include the time variable into the single-particle indices setting 1 = 
{k\, ti}. In this notation the Bethe-Salpeter equation BSE for the response function R reads: 

R(U, 23) = G(l, 3)G(4, 2) - % G(l, 5)G(6, 2)V(58, 67)i?(74, 83), (21) 

5678 

where the summation over the number indices 1, 2, . . . implies integration over the respec- 
tive time variables. The function G is the exact single-particle Green's function, and V is 
the amplitude of the effective interaction irreducible in the p/i-channel. This amplitude is 
determined as a variational derivative of the full self-energy E with respect to the exact 
single-particle Green's function: 

w> = i §|f < 22 > 

Introducing the free response R°(14, 23) = G(l, 3)G(4, 2) the Bethe-Salpeter equation (j2Tj) 
can be written in a shorthand notation as 

R = R°- iR°VR. (23) 

For the sake of simplicity, we will use this shorthand notation frequently in the following 
discussions. Since the self-energy in Eq. (J7|) has two parts £ = S+S e , the effective interaction 
V in Eq. (|2"TT) is a sum of the static RMF interaction V and time- dependent terms V e : 

V = V + V e , (24) 

where (with t\2 =t\ — £2) 

1/(14, 23) = V klk4Mk3 5(t 31 )5(t 21 )S(t u ) (25) 
is the static part of the interaction (see Eq. (!20l) ) and 

^wffrr < 26 > 

12 



contains the energy dependence. In the space of the Dirac basis ( TTTj) the amplitude V e has 
the form: 

£i(er)* /i(er) 

In order to make the Bethe-Salpeter equation (121]) more convenient for the further analysis 
we eliminate the exact Green's function G and rewrite it in terms of the mean field Green's 
function G: 

G{1, 2) = -ia kl 6 klk2 9(a kl r)e- i£ ^ T , r = h - t 2 , (28) 
where 8(t) is the Heaviside function. After a Fourier transformation in time this reads 

G klk M = ^ • (29) 

Using the connection between G and G in the Nambu form 

G- 1 (l,2) = G- 1 (l,2) + S e (l,2), (30) 
one can rewrite Eq. (12T1) as follows: 

R = R°- iR°WR, (31) 
with R°(U, 23) = (5(1, 3)G(4, 2) and W is a new interaction of the form 

W = V + W e (32) 

where 

W e (U } 23) = V e (U, 23) + zS e (l, 3)G- X (4, 2) + iG~\\, 3)S e (4, 2) - iS e (l, 3)S e (4, 2). (33) 

Since the mean field Green's function G in Eq. (|29|) is known one has a more convenient 
starting point for the solution of the Bethe-Salpeter equation ( 13~TT) than with the unknown 
exact single-particle Green's function, but the effective interaction W in Eq. ( 1321) becomes 
more complicated. The graphical representation of the Eq. (l3Ti) is shown in Fig. [2j 
In addition to the static interaction V the effective interaction W contains diagrams with 
energy dependent self-energies and an energy dependent induced interaction, where a phonon 
is exchanged between the particle and the hole. Here and hereinafter we omit the term 
iE e (3, l)E e (2, 4) in the Eq. ( 1331) by the following reasons. This term compensates the multiple 
counting of the particle-phonon coupling arising from the two previous terms in the Eq. ( 1331) . 

13 




FIG. 2: Bethe-Salpeter equation for the p/i-response function R in graphical representation. Details 
are given in Fig. Q] and the small black circle means the static part of the residual p/i-interaction 
(El. 



However, this multiple counting does not take place within time blocking approximation 
(see below) if the backward-going propagators are not taken into account. Components 
containing the backward-going propagators within lpl/i£g>phonon configurations require a 
special consideration. They have been analyzed in great detail in the Refs. 3J, |35|, [36] . In 
the present work they are fully neglected and therefore the term i£ e (3, l)£ e (2, 4) has also to 
be omitted. However, we have to emphasize, that we only neglect ground state correlations 
(backward-going diagrams) caused by the particle-phonon coupling. All the RPA ground 
state correlations are taken into account, because it is well known that they play a central 
role for the conservation of currents and sum rules. We consider that this is a reasonable 
approximation which is applied and discussed also in many non-relativistic models (see e.g. 



Refs. 



29 



30 



31 



32 



33 



34 



35 



36 



37 



38] and references therein). 



Considering a Fourier transformation of the Bethe-Salpeter equation fl3Tl) in time 



Rk 1 k 4 ,k 2 k 3 (^, £ ) — Crfcifca^S + u)Gk 4 k 2 ( £ ) 



de' 



z — ' / Lixi 

k 5 k e k 7 ks 

(34) 

one finds that both the solution of this equation R and its kernel W are singular with respect 
to the energy variables. Another difficulty arises because Eq. fl3TT) contains integrations over 
all time points of the intermediate states. This means that many configurations which are 
actually more complex than lpl/i£*Dphonon are contained in the exact response function. In 
the Ref. [34j a special time-projection technique was introduced to block the p/i-propagation 



through these complex intermediate states. It has been shown that for this type of response 



14 



it is possible to reduce the integral equation (}3"I]) to a relatively simple algebraic equation. 
Obviously, this method can be applied straightforwardly in our case too. 
Starting from the Bethe-Salpeter equation (I3TI) we divide the problem to find the exact 
response function into two parts. First we calculate the correlated propagator R e which 
describes the propagation under the influence of the interaction W e 

R e = R°- iR°W e R e . (35) 

It contains all the effects of particle-phonon coupling and all the singularities of the integral 
part of the initial BSE. Second, we have to solve the remaining equation for the full response 
function R 

R = R e - iR e VR. (36) 

Eq. fl36l) contains only the static effective interaction V and can be easily solved if R e is 
known. Thus, the main problem is to calculate the correlated propagator R e . 



D. Time blocking approximation (TBA) and subtraction method 

In order to calculate the correlated propagator R e we represent this quantity as an infinite 
series of graphs which contain mean field p/i-propagators alternated with single interaction 
acts. This can be expressed by the system of the following equations: 

R e = R°- iR°T e R , (37) 
r e = W e - iW e R°T e . (38) 

According to the main idea of the time blocking approximation we modify the integral part 



of Eq. ( 1381) making use of the time-projection operator in the form 341 ]: 



6(14,23) = 5 aki ^ ak2 5 klk3 6 k2k4 9(a kl t M )9(a kl t 2 3), (39) 

where 9{t) is the Heaviside function. This projection operator is introduced into the mean 
field propagator R° in the integral part of the Eq. ( |38l) to order the p/i-propagation in time 
and, thus, to separate in time the acts of the particle-phonon interaction from each other 
and to exclude configurations which are more complex than lpl/i®phonon. We replace in 
Eq. fl38|) the mean field propagator R° by 

R°(U, 23) -> JZ°(14, 23)6(14, 23) (40) 

15 



and obtain instead of Eq. 

r e (14, 23) = W e (U, 23) + - W e (16, 25)^(58, 67)6(58, 67)r e (74, 83). (41) 



i 

5678 



After a Fourier transformation in time we restrict ourselves to the response function 

-Rfclfc4,fc2fc 3 (^) 

OO 

Rk^MkM = ~ % J dhdhdtzdhSih - h)8(t 3 - Q6(U)e luJtl3 R(U, 23), (42) 

— OO 

which depends only on one energy variable u. 

The time projection by the operator (1391 leads to a separation of integrations in Eq. ( 1341) 
and we find an algebraic equation for the function fj4"2]) : 

Rk 1 k 4 ,k 2 k 3 (. LJ ) - Rk^kiMH^) + Rk^Mksi^WksksMkri^R-krkitokaiu), (43) 

k 5 k 6 k 7 k s 

or in shorthand notation 

R(w) = R°(uj) + R°(lu)W(uj)R(uj), (44) 

where 

W^ifc^fofoV^) = Vk^fakz + < ^fcifc 4 ,fc2fc 3 (^') (45) 

and 

^fcifc4,fcafcs( w ) = ^fe 1 fe 2 (^)4 1 fe 3 ^ 2 fe 4 - (46) 
i? / i Cl fc 2 (^) is the mean field propagator: 

£pfcM = , R ah {u) = , (47) 

Rh P (w) = — , Rha(w) = — , (48) 

£ph + & £ ah + k> 



£ P h = £ P —£h and $ is the particle-phonon coupling amplitude with the following components: 

®ph , ,hp'{ u} ) = ^s pp ' 



A 4 

lh"\Clh"h' 



io — e v + e h n — Qf 1 
^ if 

^ ( y lpp"lp'p" lpa"lp'a" \ 

p" r a" 

fj.fi* fi* fi 

^pp'^hh' | ip'pih'h y 

C<J — E p i + £h — Vt^ UJ — E p + Eh' — J J ' 



(49) 
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H M* 
lh"}dh"h' 



h „ uj-e a + e h » - tt>* 



_ Oh ~ l " 



c<j — e a " + e h — Qt 1 *—f u> — e v » + e h — 

a" p" ^ 

fi fi* fi* fi 

^aa'^hh 1 _|_ la'oHh'h 



„ l ■ ( 50 ) 

,u - e a > + e h - iiv uj - e a + e h > - U» J J 
Indices p, a and h traditionally denote the particle, antiparticle and hole types of the Dirac 
states, respectively There are, in principle, also non-zero amplitudes of the types <& p h',ha-, 
®ah',hp which cause transitions of particle-hole (ph) to antiparticle-hole (ah) pairs. However, 



in the present work we neglect them because as it has been investigated in Ref. [44| the effect 
of this kind of terms on the self-energy is very small whereas they require a lot of numerical 
effort. The reason is that these components as well as the component fl50l) contain large 
energy denominators. As it was already mentioned, the amplitudes of the types $> pp > t hh' and 
&hh', PP ' are also disregarded within our approximation. Therefore ground state correlations 
are taken into account only on the RPA level due to the presence of the V pp i^h'i Vhh', PP ' terms 
of the static interaction in the Eq. (H3l) . By definition, the propagator R(u) in Eq. (H3l) 
contains only configurations which are not more complex than lpl/i£*Dphonon. 
An important correction has to be done in the Eq. (1431) . Being adjusted to experimental 
data the RMF ground state contains effectively many correlations and, in particular, also 
admixtures of phonons. Therefore in our approach, as well as in other approaches beyond 
RPA where more complex configurations are taken into account explicitly, this admixtures 
would lead to double counting of correlations. In the present method, all the correlations 
entering through the admixture of phonons are taken care of by an additional interaction 
term: <3>(u;). Part of this interaction is therefore already contained in the static mean field 
interaction V. Since the parameter of the density functional and as a consequence the 
effective interaction V is adjusted to experimental ground state properties at the energy 
uj — 0, this part of the interaction $>(ou), which is already contained in V is given by $(0). In 
order to avoid double counting of correlations a subtraction procedure has been developed 



in the Ref. 37| were this part is removed. As a consequence one has to replace the effective 



interaction (1451) of the response equation (1431) by the function 5$: 

-> 5$(u) = - $(0). (51) 
The physical meaning of this subtraction is clear: the average value of the particle- vibration 
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coupling amplitude <3> at the ground state is supposed to be contained already in the residual 
interaction V, therefore we should take into account only the additional energy dependence, 
i.e. 5§(uj) = $(o;) — $(0) on top of the effective interaction V. Instead of Eq. (H3|) we finally 
solve the following response equation 

R = R° + R°[V + 5$]R. (52) 



E. Strength function and transition densities 

To describe the observed spectrum of the excited nucleus in a weak external field P, as for 
instance a dipole field, one needs to calculate the strength function: 

S(E) = ~ lim ImU P p(E + iA), (53) 

expressed through the polarizability Upp defined as 

n PP (u;) = P*R(u)P := £ P^iW^M^W (54) 

kik2k3k4 

The imaginary part A of the energy variable is introduced for convenience in order to obtain 
a more smoothed envelope of the spectrum. This parameter has the meaning of an addi- 
tional artificial width for each excitation. This width emulates effectively contributions from 
configurations which are not taken into account explicitly in our approach. 
In order to calculate the strength function it is convenient to convolute Eq. (1521) with an 
external field operator and introduce the density matrix variation Sp in the external field P: 

^Pjfclfe^) = RkxkiMkA^Pk-ikAi (55) 
kzki 

6p° klk2 (u) = J]^^. (56) 

Using Eq. ( 1321 we find that 5p(u) obeys the equation 

5p(u) = 5p°(uj) + R°(uj) (v + 5<S>(Lu)y p (uj) J (57) 
and the strength function is expressed as 

S(E) = — lim ImTr[P^6p(E + iA)]. (58) 

71 
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The Bethe-Salpeter Eq. (j521) gives us also the possibility to calculate the transition density 

PU = (0|<^» (59) 

for the excited state \v) at the energy Q u . In the vicinity of VL U the response function has a 
simple pole structure 

RU A M k M - P -fr^ • (60) 

which leads to 

PU 2 = 2% \l nS(W) Im 6pk ^ U + iA) ' (61) 

In order to determine the norm of the transition densities, it is convenient to rewrite Eq. (1521) 
in the following form: 

((R )- 1 (u)-V-5$(oj))r(uj) = 1 (62) 

Using Eq. ( 1601) and taking the derivative with respect to lu, we obtain the generalized nor- 
malization condition: 



N — -j- 



\w=n v 



P" = l, (63) 



P 

with the RPA norm 

In the limiting case of an energy-independent interaction, i.e. if one neglects $(tu), this 
reduces to the usual RPA normalization: 

E(w 2 -w 2 ) = 1 - ( 64 ) 



ph 



In the particle- vibration coupling model the quantity -/V^j$ in the Eq. f[6"B"j) is a non-positively 
definite matrix. Therefore, all the eigenvalues of the operator [1 — N (x;3>) L=c] _1 are less 
than or equal to 1, in analogy to the spectroscopic factor of a single-particle state. For this 
reason the sum 

J2{\pW ~ \Ph/) (65) 

ph 

is practically always less than 1. The reduction of the norm is caused by the spreading of 
the strength over lpl/i£g>phonon configurations. 
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III. NUMERICAL DETAILS, RESULTS, AND DISCUSSION 



A. General scheme of the calculations 

The method developed in the last section is applied to a quantitative description of isoscalar 
monopole and isovector dipole giant resonances in the even-even spherical nuclei 208 Pb and 
132 Sn. Our calculations are based on the energy functional with the non-linear parameter 
set NL3 [45)]. The scheme consists of three main parts: 

i) The Dirac equation for single nucleons together with the Klein-Gordon equations for meson 
fields are solved simultaneously in a self-consistent way to obtain the single-particle basis 
(Dirac basis). 

n 

ii) The RRPA equations [16| are solved with the static interaction V of Eq. fTSOj) to determine 
the low-lying collective vibrations (phonons). These two sets of particles (holes) and phonons 
form the multitude of lpl/i£g>phonon configurations which enter the particle-phonon coupling 
amplitude <5$(u;). 

iii) Eq. (|57|) for the density matrix variation 5p in the external field P is solved using 
this additional amplitude in the effective interaction V + 5<&(u). Eq. fl58l) finally allows 
to calculate the strength function corresponding to the operator P. It is found that the 
energy dependent term 5<&(u), which describes the change of the effective interaction due to 
the energy dependence of the particle- vibration coupling, provides a considerable enrichment 
of the calculated spectrum as compared to the pure RRPA. 

B. Choice of representation and basic approximations 

Eq. floTj) for 5p can be solved in various representations. In Dirac space its dimension is 
the number of p/i-pairs. In relativistic nuclear calculations it is often important to take 
into account the contributions of the Dirac sea and then the total number of ph and ah 
pairs entering the Eq. (I5T1) increases considerably with the nuclear mass number because 
of the high level density. As it was investigated in a series of RRPA calculations 



the completeness of the ph (ah) basis is very important for calculations of giant resonance 
characteristics as well as for current conservation and a proper treatment of symmetries, in 
particular, the dipole spurious state originating from the violation of translation symmetry 
on the mean field level. On the other side, the use of a large basis requires a considerable 
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numerical effort and, therefore it is reasonable to solve the Eq. (|57p in a different more 

appropriate representation. 

Two facts simplify our practical calculations: 

i) Due to the pole structure of the 5$ amplitude it turns out that the effects of particle- 
phonon coupling are only important in the vicinity of the Fermi surface and therefore one 
can restrict the number of p/i-configurations in this case much more than in the case of the 
static interaction V. 

ii) The effective interaction V, which is based on the exchange of mesons, contains only direct 
terms and no exchange terms. Therefore it can be written as a sum of separable interactions. 
It order to exploit these two effects we solve the response equation for a fixed value of the 
energy u in two steps. First we calculate the correlated propagator R e {uj) which describes 
the propagation under the influence of the interaction 8<&(u) 

R e = R° + R°5$R e . (66) 

It contains all the effects of particle-phonon coupling and all the singularities of the integral 
part of the initial BSE. Since 5$ is not separable Eq. ( 1661) has to be solved in Dirac space. 
This requires in principle for each value of u the inversion of a very large matrix. However 
the numerical effort can be reduced considerably by taking into account that the effects 
of the particle-phonon coupling are only important in the vicinity of the Fermi surface. 
Therefore the summation in the Eq. (1661) is performed only among the p/i-pairs with e p h < 
e W i n . Consequently, the correlated propagator R e differs from the mean field propagator R° 
only within this window. This approximation has been investigated numerically by explicit 
calculations with different window sizes as it is discussed below. It is important to emphasize 
that the number of ph- and a/i-configurations taken into account on the RRPA level has 
to be considerably larger in order to obtain convergence in the centroid positions of giant 
resonances as well as to find the dipole spurious state close to zero energy. 
In the second step, we have to solve the remaining equation for the full response function 
R(u): 

R = R e + R e VR. (67) 

In contrast to the Eq. fl52|) this equation contains only the static effective interaction. 
Now we exploit the fact that the one-boson exchange (OBE) interaction discussed in Ap- 
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pendix O is separable in momentum space. It can be expressed as 

Vfcifc4,fc 2 fc3 = dc Qklk 2 QkTki ( 68 ) 

c 

where the channel index c = (g, m) is given by the momentum q transferred in the exchange 
process of the corresponding meson labeled by the index m. The parameters d c are given by 
the meson propagators. We now can use the well known techniques of the response formalism 
with separable interactions (see, for instance, Ref. ^M\)- We define 

R cd (u) = ^ Qkik2 R kxkAMh 3 {^)Qk'X ( 69 ) 

k^k^ki 

and find from Eq. (loTD 

R cd = R lc> + R cc" d c" R c"c' = R lc> + (R e dR) cc >. (70) 

c" 

This equation can be solved by matrix inversion 

R = (l - R e d^j R e . (71) 

Finally we transform Eq. (137)1 to the momentum space, which is equivalent to using the 
external field as an additional channel c = p, and obtain from Eq. (I6T?|) the full response 

5p c = 5p e c + Y,R e cc'd c ,5p c > (72) 

d 

and the polarizability (J51|) : 

U PP = P^RP = P^R e P + J2 6 Pc d c S Pc = R £ PP + E^X(l - Red )jp e c>- (73) 

c cc' 

The rank of vectors and matrices entering the Eq. (1721) is determined by the number of 
mesh-points in g-space and the number of meson channels. In particular it does not depend 
on the mass number of the nucleus. In realistic calculations the rank of the Eq. (1711) in 
the momentum space is around 500 and, obviously, the numerical effort does not depend 
considerably on the total dimension of ph and ah subspaces. However, if one stays in the 
Dirac basis the latter dimension is exactly the rank of arrays in the Eq. (1431) . One should 
bear in mind that in practice both subspaces (ph and ah) are truncated at some energy 
differences e p h and e a h which are large enough so that a further increase of these values 
does not influence the results. In light nuclei with relatively small level density the total 
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dimension of ph and ah subspaces is comparable and could be even smaller than the rank of 
the Eq. (1721) and, therefore, working in the Dirac basis is more preferable. For heavy nuclei 
the dimension of the Dirac basis becomes huge and, therefore, using of the momentum space 
is more justified. 



C. Numerical details 



To ensure numerical properness of our codes the response equation for dp has been solved 
both in momentum space f!72|) and in Dirac space fl57j) and identical results have been 
obtained. Both Fermi and Dirac subspaces were truncated at energies far away from the 
Fermi surface: in the present work as well as in the Ref. 44( we fix the limits e p h < 100 MeV 
and e a h > —1800 MeV with respect to the positive continuum. A small artificial width of 200 
keV was introduced as an imaginary part of the energy variable u to have a smooth envelope 
of the calculated curves. The energies and amplitudes of the most collective phonon modes 
with spin and parity 2 + , 3~, 4 + , 5~, 6 + have been calculated with the same restrictions 
and selected using the same criterion as in the Ref. 44] and in many other non-relativistic 
investigations in this context. Only the phonons with energies below the neutron separation 
energy enter the phonon space since the contributions of the higher-lying modes are found to 



be small. Test calculations in the framework of the approach |37l . |38| without the restriction 
of the phonon space by the energy resulted in the small deviation of the strength function as 
well as the change of the mean energies and widths of the resonances comparable with the 
smearing parameter (imaginary part of the energy variable) used in the calculations. This 
is the natural result because the physical sense of this parameter is to emulate contributions 
of the configurations which are not taken into account explicitly. 

Because of the pole structure of the particle-phonon coupling amplitude (l4"9j) and (!50|) its 
contributions to the final result decrease considerably when we go away from the Fermi 
surface. Therefore this coupling has been taken into account only within the ph-energj 
window e p h < 30 MeV around the Fermi surface. It has been checked that further increase 
of this window does not influence considerably the results. 

It is important to note that although a large number of configurations of the lpl/i®phonon 
type are taken into account explicitly in our approach, nevertheless we stay in the same 
ph (ah) space as in the RRPA, therefore the problem of completeness of the phonon basis 
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does not arise and, so, the phonon subspace can be truncated in the above mentioned way. 
Another essential point is, that on all three stages of our calculations the same relativistic 
nucleon-nucleon interaction V (|20|) has been employed. The vertices 7fc lfca (fl9|) entering the 
term 5Q(uj) are calculated with the same force. Therefore no further parameters are needed. 
Our calculational scheme is fully consistent. 



The subtraction procedure (T5TT) developed in the Ref. 37J for the self-consistent scheme has 
been incorporated in our approach. As it was mentioned above, this procedure removes the 
static contribution of the particle-phonon coupling from the p/i-interaction. It takes into 
account only the additional energy dependence introduced by the dynamics of the system. 



It has been found in the present calculations as well as in the calculations of the Ref. 38] that 
within the relatively large energy interval (0 - 30 MeV) the subtraction procedure provides 
a rather small increase of the mean energy of the giant dipole resonance (0.8 MeV for lead 
region) and gives rise to the change by a few percents in the sum rule. This procedure restores 
the response at zero energy and therefore it does not disturb the symmetry properties of 
the RRPA calculations. The zero energy modes connected with the spontaneous symmetry 
breaking in the mean field solutions, as for instance the translational mode in the dipole case, 
remain at exactly the same position after the inclusion of the particle-vibration coupling. 
In practice, however, because of the limited number of oscillator shells in our calculations 
this state is found already in RRPA without particle-vibration coupling at a few hundreds 
keV above zero. In cases, where the results depend strongly on a proper separation of this 
spurious state, as for instance for investigations of the pygmy dipole resonance in neutron 
rich systems, we have to include a large number of ph -configurations in the RRPA solution. 

D. Isoscalar monopole and isovector dipole resonances in 208 Pb and 132 Sn 

The calculated strength functions for the isoscalar monopole resonance in 208 Pb and 132 Sn 
are shown in Fig. [3] The fragmentation of the resonance caused by the particle-phonon 
coupling is clearly demonstrated although the spreading width of the monopole resonance 
is not large because of a strong cancellation between the self-energy diagrams and diagrams 
with the phonon exchange (see Fig. [2]). This fact has also been discussed in detail in Refs. 
(3, [30] and it is not disturbed by the subtraction procedure floTl) because this cancellation 
takes place as well in $>(u) as in $(0). 
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FIG. 3: Isoscalar monopole resonance in 208 Pb and 132 Sn obtained within two approaches: RRPA 
(dashed line) and RRPA with particle-phonon coupling RRPA-PC (solid line). Both computations 
have been performed with relativistic Hartree (RH) mean field and employ NL3 parameter set for 
RMF forces. 

In order to compare the spreading of the theoretical strength distributions with experimental 
data we deduce mean energies (E) and widths parameters T by fitting our theoretical strength 
distribution in a certain energy interval to a Lorentz curve in the same way as it has been done 
in the experimental investigations. The mean energies and the widths parameters obtained 
in this way are shown in the Table [B As experimental data we display the numbers adopted 



in the Ref. 48] from the evaluation of a series of data obtained in different experiments for 
the isoscalar monopole resonance in 208 Pb. 



25 



TABLE I: Lorentz fit parameters of isoscalar EO strength function in 208 Pb and 132 Sn calculated 
within RRPA and RRPA extended by the particle-phonon coupling model (RRPA-PC) as compared 
to experimental data. The fit has been carried out in the interval from B n to roughly 20 MeV. 





<E> (MeV) 


r (MeV) 


RRPA 


14.16 


1.71 


208 Pb RRPA-PC 


14.05 


2.36 


Exp. [48] 


13.73(20) 


2.58(20) 


RRPA 


16.10 


2.63 


132 Sn RRPA-PC 


16.01 


3.09 



The calculated photoabsorption cross sections 

167r 3 e 2 

oei{E) = -^—E Sei(E) (74) 

for the isovector dipole resonance in 208 Pb and 132 Sn are given in the Figs. H] and [5] re- 
spectively. The left panels show the results obtained within the non-relativistic approach 
with Woods-Saxon (WS) single-particle input and Landau-Migdal (LM) forces described in 



Ref. 



381 ] . They are compared in the right panel with the relativistic fully consistent th eory 



developed in the present work. The underlying Lagrangian has the parameter set NL3 45]. 
RPA calculations are shown by the dashed curves, RPA-PC calculations - by the thick solid 
curves. In Fig. H]we have displayed experimental data with error bars taken from Ref. 501 ] . 
To make the comparison reasonable calculations within the non-relativistic framework have 
been performed with box boundary conditions for the Schrodinger equation in r-space which 
ensures completeness of the single-particle basis. 

The corresponding Lorentz fit parameters in the two energy intervals: B n —25 MeV and 0—30 
MeV (B n is the neutron separation energy) are included in Table [III and they are compared 



with the data of Ref. 



49 



501 ] . We notice that the inclusion of particle-phonon coupling 



in the RRPA calculation induces a pronounced fragmentation of the photoabsorption cross 
sections, and brings the width of the GDR in much better agreement with the data, both 
for 208 Pb and 132 Sn. 
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FIG. 4: The isovector El resonance in 208 Pb. Details are given in the text. 

The fragmentation of the resonance introduced by the particle-phonon coupling is clearly 
demonstrated in both cases. Also, one finds more or less the same level of agreement be- 
tween theory and experimental data for these two calculations. In the case of the isovector 
El resonance in 132 Sn this is, however, not so clear because the cross section and the inte- 



gral characteristics of the resonance obtained in the experiment of Ref. 49j are given with 
relatively large error bars. In 208 Pb our self-consistent relativistic approach reproduces the 
shape of the giant dipole resonance much better than the non-relativistic one although the 
whole resonance is about 0.5 MeV shifted to lower energies with respect to the experiment. 
As one can see from the Fig. H] and Table HIl we observe some shift already in the RRPA 



calculation, which is determined by the properties of the NL3 
forces, for instance, the use of the density dependent versions 



orces. Improvement of the 



51 



52| of the RMF should 
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FIG. 5: The same as in Fig. [3] but for 132 Sn. 
bring the El mean energy in better agreement with the data. 

However, there is an essential difference between the fully self-consistent relativistic calcu- 
lations of the present work and the non-relativistic approach: in non-relativistic approach 
discussed in Ref. [38J one introduces on all three stages of the calculation phenomenological 
parameters, which have to be adjusted to experimental data: first, the Woods-Saxon param- 
eters as, for instance, the well depth are varied to obtain single-particle levels close to the 
experimental values, second, one of the parameters of the Landau-Migdal force is adjusted 
to get phonon energies at the experimental positions (for each mode) and, third, another 
Landau-Migdal force parameter is varied to reproduce the centroid of the giant resonance. 
Although the varying of the parameters is performed in relatively narrow limits, it is nec- 
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TABLE II: Lorentz fit parameters in the two energy intervals: B n — 25 MeV and — 30 MeV, for 
the El photoabsorption cross sections in 208 Pb and 132 Sn, calculated with the RRPA, and with the 
RRPA extended to include the particle-phonon coupling (RRPA-PC), compared to data. 



B n - 25 MeV - 30 MeV 







<E> 


r 


EWSR 


<E> 


r 


EWSR 






(MeV) 


(MeV) 


(%) 


(MeV) 


(MeV) 


(%) 




RRPA 


13.1 


2.4 


121 


12.9 


2.0 


128 


208pk 


RRPA-PC 


12.9 


4.3 


119 


13.2 


3.0 


128 




Exp. [50] 


13.4 


4.1 


117 






125(8) 




RRPA 


14.7 


3.3 


116 


14.5 


2.6 


126 


132 Sn 


RRPA-PC 


14.4 


4.0 


112 


14.6 


3.2 


126 




Exp. [49] 


16.1(7) 


4.7(2.1) 


125(32) 









essary to obtain realistic results. In contrast, in the relativistic fully consistent approach 
developed in the present work no adjustment of additional parameters is made. Of course, 
the underlying energy functional has been determined in a phenomenological way by a fit 
to experimental ground state properties of characteristic nuclei. However, it is of universal 
nature and the same parameters are used for investigations of many nuclear properties all 



over the periodic table. The predictive power of this scheme is therefore much higher_ 
that of the present semi-phenomenological approach discussed, for instance, in Ref. 



;han 



The effect of the particle-vibration coupling on the low-lying dipole strength around the 
neutron threshold within the presented approach can be also seen in the Figs. I4"f5l But this 
phenomenon requires the detailed investigation and is considered in the special publication 



53]. 
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IV. SUMMARY AND OUTLOOK 



In order to describe the damping of collective excitations in nuclei covariant density func- 
tional theory is extended by particle-vibration coupling in a consistent way. Starting from 
a relativistic energy functional E[p] the self-consistent RMF equations are solved and the 
static part of the nucleon self-energy £ = SE/Sp is found. In a second step the low-lying 
collective phonons are determined in the framework of relativistic RPA with the effective 
p/i-interaction V = 5T*/5p. Using particle-phonon vertices deduced from the same interac- 
tion the self-energy £ is extended by an energy dependent part S e (e). The Bethe-Salpeter 
equation derived from the full self-energy S = £ + S e is formulated in the ph-basis of the 
Dirac eigenstates. Using the time blocking approximation, which allows the truncation to 
lpl/iCg>phonon configurations and a subtraction procedure, which avoids double counting of 
correlations it is shown that the resulting effective interaction V + 5§(uj) contains besides 
the static part V an energy dependent correction term 5Q(u), which takes into account the 
coupling of the particles to collective phonons. 

This method is applied to the computation of spectroscopic characteristics of nuclear excited 
states in a wide energy range up to 30 MeV for spherical nuclei with closed shells. An 
equation for the density matrix variation is formulated and solved in the Dirac space as well 
in the momentum space. The particle-phonon coupling amplitudes of collective vibrational 
modes below the neutron separation energy have been calculated within the self-consistent 
RRPA using the parameter set NL3 for the Lagrangian. The same force has been used for 
the static part of the effective p/i-interaction V and for the evaluation of the particle-phonon 
vertices 7^,. Therefore a fully consistent description of giant resonances is performed. 
Noticeable fragmentation of the isoscalar monopole and isovector dipole giant resonances in 
208 Pb and 132 Sn is obtained due to the particle- vibration coupling. This leads to a significant 
spreading width as compared to simple RRPA calculations. This is in agreement with 



experimental data as well as with the resu 



phenomenological approaches of Refs. 
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38]. 



ts obtained within the non-relativistic semi- 



Thus in the present work a description of nuclear many-body dynamics including complex 
configurations is realized within an approach which is (i) fully consistent, (ii) based on 
relativistic dynamics, (iii) universally valid for nuclei all over the periodic table, and (iv) 
based on a modern covariant density functional, which has been applied with great success 



30 



for many nuclear properties all over the periodic table. 

So far, this extended density functional theory for complex configurations has been formu- 
lated only for closed shell nuclei. To expand the field of applications to nuclei with opened 
shells it is, of course, necessary to include pairing correlations into the approach discussed in 



the present paper. This can be done in the manner of Ref. [37] where pairing correlations and 
particle-vibration coupling are taken into account on an equal footing in a non-relativistic 
framework using techniques of a generalized Green's function formalism. Work on this di- 
rection is in progress. 
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APPENDIX A: NON-LINEAR MESON POTENTIALS 

So far we started from a Lagrangian ([T]) containing the nucleons as Dirac particles and mesons 



providing an interaction between them 



46, 



4=71 ] . as it is the starting point of most of the 



relativistic investigations in nuclear physics. It has been recognized, however, very early 54] 



that this model does not provide an adequate description of the surface properties of realistic 
nuclei. The incompressibility is too high and the deformations are too small Therefore 
over the years the relativistic models have been considerably improved by introducing an 
effective density dependence, as it is in accordance with the concept of density functional 
theory. Boguta and Bodmer proposed already in 1977 to introduce a non-linear coupling 
between the q-mesons replacing the mass term ^m 2 a a 2 in the Lagrangian ([I]) by a non-linear 



meson potential 54] 



U{a) = -m a a + —a + — a . (Al) 
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This procedure was used because it did not destroy the renormalizability of the model. In 
these early days relativistic models of the nucleus were considered as fully fledged quantum 
field theories 47] and therefore non-linear meson couplings provided an elegant extension of 
the model containing additional parameters introducing an effective density dependence. 
In fact, models with non-linear meson couplings have been used very extensively in the liter- 
ature. Very successful parameterizations have been developed on this basis as, for instance, 
NL1 [55] and NL3 45], which contain only non-linear couplings between the cr-mesons. Later 



on one also considered non-linear couplings between the w-mesons |56|,|57|, however without 
improving the agreement between theoretical results and experimental data in finite nuclei. 
In the present investigations we consider the parameter set NL3, which has been used with 
great success for the description of symmetric nuclear matter and finite nuclei with closed 
shells deformed nuclei, rotating nuclei [ijj and for relativistic RPA calculations of giant 
resonances and low-lying states [58]. In the most successful applications pairing correlations 
are taken into account in the framework of relativistic Hartree-Bogoliubov (RHB) theory 



with Gogny type interactions in the pp-channel 60]. For a recent review see Ref. |6j. 
The effective interaction entering the static part of the nucleon self-energy is in this case 
given by meson exchange potentials. Without the non-linear couplings these interactions are 
of Yukawa type with a range determined by the mass of the corresponding mesons and with 
a relativistic structure determined by the quantum numbers of the mesons. In the simplest 
case of the crtu-model we have 

V(l, 2) = -glD a (ru r 2 )/3«/3^ + r 2 )(l^ - a ® a V>) (A2) 

with the Dirac matrices (3 and a and with the propagator in r-space 

]_ e -m m \r 1 -r 2 \ 

D m (r 1 ,r 2 ) = - : r , for m = a,u. (A3) 

Aix |r*i — r 2 \ 

More realistic applications contain in addition a p-meson and the photon. These potentials 
are obtained from the Klein-Gordon equations by elimination of the meson degrees of freedom 
neglecting retardation. 

In the non-linear case the static Klein-Gordon equation for the a-meson has the form 



- Aa + U'(a) = -g a (W) , (A4) 

with 

U'(a)=mla + g 2 a 2 + g 3 a 3 . (A5) 
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This is a nonlinear equation, which does not allow for an analytic solution. The numerical 
solution of Eq. (1A4I) cr(r) provides us with the mean field of scalar type 

Z.(r)=g a(r)- (A6) 



Re 



ativistic RPA is obtained as the small amplitude limit of time-dependent Hartree theory 



161 ] . In this case the relativistic single-particle density matrix pkk> is expanded around the 



static solution po- The linearization of the Klein-Gordon equation (IA4I) leads to 

(-A + M(r))5a = -gJp s , (A7) 

where 5p s is the scalar density, with M(r) depending on the coordinate r 

M{r) = U"(a (r)) = m 2 a + 2g 2 a (r) + 3g 3 a 2 (r). (A8) 

where c (r) is the cr-field of the static solutions. In momentum space this leads to a non-local 
propagator, which is the solution of the integral equation 

q 2 D a (q,q') + J^M(q- q")D a (q", q') = (2n) 3 5(q - q'), (A9) 

where M(q) is the Fourier transform of M(r): 

M{q) = Jd 3 re- iqr M(r). (A10) 

The a part of the effective ph-interaction used in the RPA calculations has therefore the 
form 

Vf(l, 2) = -gl^^ J d ^0 r e^-^DM, q') (AH) 

and similar terms for the other mesons and the photon. We use this interaction in the 
solution of the relativistic RPA equation for the calculation of the collective phonons, as well 
as in Eq. f|T9|) for the determination of the corresponding particle-phonon vertices 7^ and 
in Eq. (|36|) for the response function i?(14, 23). 
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APPENDIX B: DENSITY DEPENDENT MESON VERTICES 



In recent years it has been recognized that the relativistic models are by no means fully 
fledged quantum field theories. Instead of that they are effective field theories, which provide 
the basis of covariant density functional theory. Renormalizability is not important. In the 
spirit of density functional theory it is reasonable to abandon the non-linear couplings and to 
introduce instead of that coupling parameters g a (p(r)), g w {p{r)) and g p (p(r)), which depend 



on the baryon density 
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62 



63 



64|. 



In this case the static self-energy contains rearrangement terms, i.e. terms depending on the 
derivative of the coupling constant with respect to the density. We find 



E a (r) = g a a{r) 

E (r) = gjjJoir) + g P T 3 p (r) + V c (r) + E K (r) 



(Bl) 
(B2) 



with the rearrangement term 



£» = g' a (p(r)) Ps (r)a(r) + g'Mr))p{r)u (r) + g' p (p(r)) p T {r) p {r) , (B3) 

where p = p p + p n and pr = P P — Pn are the isoscalar and isovector baryon densities and 
where po is the time-like component of the p-field. Here we have used the fact, that we have 
time-reversal invariance in the ground state and that the currents vanish in this case. 
The effective ph-interaction is obtained as the derivative of the self-energy with respect to 
the density operator. We find for the a-exchange 

V*(l, 2) = -^^g a {l)g a {2)D a {r u r 2 ) - J^n^n - r 2 ) 

- {^hWg c {l)Ut) + 1 (1) /3 (2) ^(1)^(2) + l«l( 2 )^ (T (l)^(2)} J D (T (r 1 ,r 2 ) (B4) 



with 



Un) = {(/3 (1) K 2 ) + i(V 2) K(i) + i (1) i (2) ^(i)P,(n)} x 

x / d 3 r ^(p(r))L) CT (ri,r)p s (r) 



(B5) 



and similar terms for the lu- and p-exchange. Here we have used the abbreviations: 

g a (l) := g a {p{ri)) 
9aiX) ■= 9 ! a(p( r i))ps{ri)- 



(B6) 
(B7) 
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APPENDIX C: SOLUTION IN MOMENTUM SPACE 



First we concentrate on the static interaction V of Eq. (jA2j) neglecting for the moment the 
non-linear meson coupling. We can transform the one-boson exchange (OBE) interaction 



(1A2I) to momentum space and neglect retardation: 



v(i, 2) = ±Yj& f 70 I e^- r ^r-(i)^( g )r-(2) (CI) 

where the -l-sign holds for the vector mesons and the —sign for the scalar mesons. The index 
m denotes the set of mesons m = {a, u, p, e} - the isoscalar-scalar cr-meson, the isoscalar- 
vector a;-meson, the isovector-vector p-meson and the photon, which carry the respective 
components of the interaction. Summation is implied also over the repeated Greek indices. 
Dm(q) is a boson propagator which is usually taken in the RMF calculations in the simplified 
Yukawa form neglecting formfactors and retardation effects: 

D m (q) = - 1 2 , m = a, (C2) 

D%{q)= " 2 , m = u,p,e. (C3) 

Introducing the channel index c = (q, m), which combines the meson index m with the 
momentum transfer q, we can express V as a sum of separable terms 

V{l,2) = Y,dcQ {c) H)Q {c)] {2) (C4) 



with 



and the propagator 



Q^ = Q( q ,m)=T™e^ (C5) 



d c = d(c l ,m) = ±g 2 m D%( q ). (C6) 
The matrix elements of V in Dirac space have the form 

C 

= ±Y,9 2 m I ^(^|r^ r |A: 2 )^( ? )(A: 4 |r-e-^|A:3). (C7) 

For non-linear meson-couplings this expression has to be somewhat modified, because the 
propagator in momentum space (IA110 is non-diagonal in this case. We have D^(q, q') and 
the values d c have to be replaced by matrices d cc > in this case. 

35 



APPENDIX D: RESPONSE FORMALISM IN SPHERICAL NUCLEI 



In the spherical case angular momentum coupling reduces the numerical effort considerably. 
The response equation for Sp reads in this case: 



+ Mfkiktteka) ( U ) 6 P'kk 4 ) (<*>)■ (Dl) 

(^3^4) 

It contains the matrix elements of the static effective interaction V J and 5Q j (uj) is the 
change of the effective interaction due to the energy dependence of the particle-vibration 
coupling. 

For the meson exchange potentials we obtain after coupling to good angular momentum J: 



(feife4,fe 2 fe 3 ) 2 J+ 1 ^ K ' / (2tt) 3 

LSm J ' 







x -^^(ks\\h(qr)[T^Y L } J \\h), (D2) 



in 



where index S = (0, 1) denotes the spherical component of the Pauli matrix entering the 
vertices The factor (— ) s indicates that the space-like parts of the interaction (current- 
current interactions) have the opposite sign as the time-like part. The matrix elements ( 1D2I ) 
are a sum of separable terms. The dimension of the matrices to be inverted is given by the 
number of separable terms. 

The interaction V is treated in the response equation Eq. fl72l) . After a transformation of 
this equation to momentum space we find 

00 

(47r) ^-^ / i\S' f 9 dq r e \j mm > . , g , j , , 

± 2~TTT ^ ( ' / l^W L ^ L ' S ' ^ q '^ \^ i —2 6p L , sl {q;u), (D3) 

L'S'm' J \ ) H 



where 



with 



4t^' E ^i5( fclfc2 )(^)< ) i 4 , fc2fc 3)(^)^ {fc 3 fc4 )(^) (D4) 

(fcifc 2 fc3fc4) 



QiScwta)fa) = (^lljL(gr)[r-F L ] J ||fc 2 ) (D5) 
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is a Fourier transform of the correlated propagator R^ J which is determined in the Dirac 
basis by the angular momentum coupled version of Eq. (1661) : 

^(fct ) fc a fcs)H = R(k lk2 )(uj)5 {klk3) 5 {k2k4) +R (klk2) (u) <y$ ffcifc8,fc2fc B )( a; )- R (^fc4 ) fc 6 *3)( a; )- ( D6 ) 

(k 5 k G ) 

The free term on the right side of the Eq. ( 1D3I) has the following form: 

(kik2k3k4) 



where 



P( klk ,) = {ki\\P J \\k 2 ). (D8) 



The spectrum of nuclear excitations in the external field P with the multipolarity J is 
therefore determined by the strength function: 

S J (E) = -- lim Im U J (E + iA), (D9) 

expressed through the following polarizability U j (uj): 

oo 

n>) = n<^M±^ 4j£«^ teu )^^^ (r ,„), (dio) 







where 

Thus, one can see that it is convenient to solve our problem as well as the RRPA problem 
in the momentum space since the rank of vectors and matrices entering the Eq. (ID 31) is 
determined by the number of mesh points in g-space and the number of meson channels. In 
particular it does not depend on the mass number of the nucleus. In realistic calculations 
the rank of the Eq. ( 1D3I) in the momentum space is around 500 and, obviously, the numerical 
effort does not depend considerably on the total dimension of ph and ah subspaces. However, 
if one stays in the Dirac basis the latter dimension is exactly the rank of arrays in the 
Eq. fIDip . One should keep in mind that in practice both subspaces are truncated at some 
energy differences e p h and e a h which are large enough so that a further increase of these 
values does not influence the results. In light nuclei with relatively small level density the 
total dimension of ph and ah subspaces is comparable and could be even smaller than the 
rank of the Eq. ( 1D3I) and, therefore, working in the Dirac basis is more preferable. For heavy 
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nuclei the Dirac basis becomes huge and, therefore, using of the momentum space is more 
justified. 

However, even if one considers the Eq. (1D3j) in the momentum representation, Eq. (1D6I) for 
the correlated propagator, nevertheless, has to be solved in the Dirac basis. This equation 
contains the particle-phonon coupling amplitude $ J with the following matrix elements 
coupled to angular momentum J: 

H Jn h" p > n 

r fl p* [1 [A* 

°x p ,xp r-s—y 7(p P '')7( p y) v-^ 7(p Q '')7(p' a ") 



, $ ""p'Xp fsr^ >(pp") Hp'p") 

(W) 2j p + l\<^f uj - e p „ + e h - Q» u-e all +e h -^ 



p 



+ (-l) J+J » 



Jp Jh -J I A x ^ l(pp')l(hh') \ L > '(p'p)'(h'h) 

3h> Jp' \ ^ ^ - £ p' + e h ~ ^ uj - s p + e h > - 



(D12) 

where x fc denotes the relativistic quantum number set: x fc = (2j k + l)(/ fc — and Jn. lk2 \ = 
(^i||7 M ll^2) denotes the reduced matrix element of the particle-phonon coupling amplitude. 
The backwards going components are found through the symmetry relations: 

^aW") = (- i y h+jp+jh ' +jp ^u,^-^- ( Dis ) 

But this problem is too expensive numerically to be solved in the full Dirac basis. Due to the 
pole structure of the $ amplitude it is naturally to suggest that particle-phonon coupling 
effect is not important quantitatively far from the Fermi surface. So, in the numerical 
calculations an energy window e W i n was implied around the Fermi surface with respect to 
ph energy differences e p h so that the summation in the Eq. (1D6j) is performed only among 
the ph pairs with e ph < £ W in- Consequently, the correlated propagator differs from the 
mean field propagator only within this window. This approximation has been investigated 
numerically by direct calculations with different values of e W i n as it is discussed in Sec. IIHI 
It is important to emphasize that many ph and ah configurations outside of the window are 
taken into account on the RRPA level that is necessary in order to obtain the reasonable 
centroid positions of giant resonances as well as to find the dipole spurious state close to 
zero energy. 
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